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Abstract 

We present a high-statistics Monte Carlo determination of the exponent 7 
for self-avoiding walks on a Manhattan lattice in two dimensions. A conser- 
vative estimate is 7 ^ 1.3425 ± 0.0003, in agreement with the universal value 
43/32 on regular lattices, but in conflict with predictions from conformal field 
theory and with a recent estimate from exact enumerations. We find strong 
corrections to scaling that seem to indicate the presence of a non-analytic ex- 
ponent A < 1. If we assume A = 11/16 we find 7 = 1.3436 ± 0.0003, where the 
error is purely statistical. 



1 Introduction 



The self-avoiding walk (SAW) is a model which describes the universal properties of 
flexible chain polymers in a good solvent in the dilute regime. A simple but intriguing 
modification has been recently introduced to study polymers with an intrinsic orien- 
tation [Q^]. This orientation could be due to the presence of dipole moments on the 
monomers of the chain or to an ordering in the sequence of monomer constituents. 
On the lattice one considers SAWs with a short-range interaction between different 
steps of the walk according to their relative orientation. The partition function is 
simply 

z N = Y,e Pama+0pmp , (1.1) 
M 

where m p and m a are respectively the number of parallel and antiparallel interactions 
and the sum extends over all SAWs of length N. 

This model has a rich phase diagram |3|-|7] and can be analyzed theoretically [|T],0 
by mapping it into a complex 0(n) model in the limit n — > 0. In two dimensions, the 
theory with a repulsive interaction between parallel bonds, i.e. with (3 a = 0, (3 P < 0, 
was analyzed using conformal invariance techniques 0. It was shown that the new 
interaction is truly marginal, giving rise to a line of fixed points. The main conse- 
quence is that the partition-function exponent 7 should vary continuously with the 
strength of the orientation-dependent interaction. The exponent 7 is defined from the 
asymptotic behaviour of the partition function zn, which should scale asymptotically 
as 

ZN = A^N^ 1 . (1.2) 

Here A and // are non-universal constants, while 7 is an exponent that, in the absence 
of the orient at ion- dependent interaction, is expected to be universal: it should not 
depend on the details of the interaction and it should assume the same value for 
any two-dimensional regular lattice. Using Coulomb-gas techniques, Nienhuis || 
predicted 7 reg = 43/32, a value that has been confirmed to high-precision by many 
numerical computations, see e.g. Ref. 0. On the other hand, for interacting oriented 
SAWs with (3 = 0, (3 P < 0, 7 should be a function of (3 P . Unfortunately, conformal 
field theory does not provide definite numerical estimates, although it predicts that 
7(/3 p ) should decrease monotonically as (3 P — > —00. 

The square-lattice model was studied by exact enumerations in Ref. []3| and by 
transfer-matrix techniques in Ref. Jl0| finding a very tiny dependence of 7 on (5, 



thereby supporting the field-theory analysis. When parallel interactions are forbidden, 
i.e. for P p = —00, Ref. || finds using unbiased approximants, 

0.006 < 7 re g- 7 < 0.013 . (1.3) 

The systematic uncertainty due to the extrapolations of two different series is taken 
into account in this range. The evidence for a non-universal behaviour is not over- 
whelming, keeping also into account that several problems affected the results of the 
analysis. As the authors report, they found small shifts of the critical fugacity /i c from 
the value it assumes for the ordinary SAW, in contrast with the theoretical result that 
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He should not depend on (3 P . More worryingly, an analysis using biased differential 
approximants with fixed critical fugacity gave a smaller prediction for 7 reg — 7, al- 
though with much less confidence since most of the approximants were defective. The 
analysis of Koo | I0| , based on strips of width < 8, was less precise. For the largest 



value of (3 P that was analyzed, (3 P = —3, he finds 7 reg — 7 = 0.018 ±0.012, a difference 
that is barely significant. 

The evidence provided by Refs. |3|,[IU| for a non-universal behaviour of 7 was not 



conclusive and this spurred many workers to improve the result and/or to investigate 
the problem using different methods. The transfer-matrix analysis was improved in 
Ref . || : using larger strips they did not find any evidence of a non-universal behaviour 
and interpreted previous results as due to short-series effects and to the small size of 
the strips. Other field-theoretical predictions were also tested. Refs. P,|TT|] computed 



by Monte Carlo and exact-enumeration methods the mean value of m p for ordinary 
SAWs: they find that (m p ) converges to a constant as iV — > 00, in contrast with 
the field-theoretical prediction (m p ) ~ logiV. Recently the behaviour of (m p ) on a 
cylinder was determined by a Monte Carlo simulation ||12|| , finding also in this case a 
result in disagreement with the field-theory predictions. 

Recently it was shown [13"] that Cardy's original argument implies that the ex- 



ponent 7 for the Manhattan lattice should also be different from the exponent 7 reg . 
Indeed on this lattice a SAW is oriented by default and parallel interactions are auto- 
matically suppressed. From the analysis of long exact-enumeration series the authors 
of Ref. |L3j report 

7re g - 7 = 0.0053 ± 0.0030 . (1.4) 

The effect is extremely small. It differs from zero by less than two error bars. The 
evidence for 7 7^ 7 re9 is therefore not overwhelmingly persuasive, and the theoretical 
importance of the problem asks for further investigations. Indeed one can suspect 
that the small deviation is simply a systematic effect due to the corrections to scaling 
that are not completely taken into account by the analysis. 

We have therefore decided to investigate the problem by means of a Monte Carlo 
simulation, computing the exponent 7 for SAWs on a Manhattan lattice. The ad- 
vantage is that we are able to work with very long walks (N < iV max = 32000) and 
therefore to reduce the unknown systematic uncertainty due to the extrapolation 
iV -> 00. 

Our simulations were performed with two different algorithms. The first one, the 
join-and-cut algorithm is a dynamic Monte Carlo algorithm that works in the 



ensemble of couples of walks with fixed total length. The algorithm is at present 
the best one to compute the exponent 7 since the autocorrelation time in CPU units 
scales as iV L6 , while for other algorithms it behaves no better than N~ 2 . The second 
algorithm is a variant of the pruned-enriched Rosenbluth method (PERM) ||15|| . This 
is a growth algorithm. Asymptotically, it is slower than the join-and-cut algorithm 
in generating independent configurations. The computer time to generate an inde- 
pendent configuration scales as N 2 . But the constant in front of N 2 is very small, for 
instance the number of monomer additions needed to obtain one independent config- 
uration was numerically found to increase as ~ 0.008iV 2 for N > 10000. Therefore 
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it is possible that the PERM is the most efficient one even for quite long walks. As 
we shall discuss, with a clever improvement, the Markovian anticipation, the PERM 
is more efficient in providing estimates of 7 than the join-and-cut algorithm as long 
as the length of the sampled walks is less than 10 4 . Only for longer walks the join- 
and-cut algorithm is faster. Another advantage in the present context is that it gives 
directly, together with the estimate of the partition sum for chains of length iV max 
and without extra cost, also all partition sums for shorter chains. These estimates for 
different N are not independent, but just because of this fact they are particularly 
useful for estimating 7. 

The main sources of systematic errors in our analysis are the corrections to scaling, 
assumed to be of the form 

Z " -1 + ^ + ^7 + - (1-5) 



A^ N y-l N A N Ai 

with A < Ai < • • -. All numerical evidence for SAWs on regular lattices (square, 
honeycomb and triangular) indicates that the leading correction is the analytic one 
|T6Hl8|. On the other hand Saleur 0] predicted A = 11/16, a result that was 
confirmed in numerical work on lattice trails [^,|2T|].[| Why this exponent does not 
show up for SAWs on regular lattices is completely unclear. 

For the Manhattan lattice, Ref. |13| showed that the enumeration data are very 



well fitted using an Ansatz with no non-analytic terms[| with A < 1. On the other 
hand a naive fit of our Monte Carlo data would indicate just the opposite: the results 
are well fitted assuming a correction-to-scaling exponent of order A « 0.5 — 0.7. Of 
course one should not take this indication too seriously — our data are not precise 
enough for a serious attempt to determine A — but it is fair to say that A = 11/16 
is our preferred value. If we assume A = 11/16 we obtain 

7 = 1.3436 ±0.0003, (1.6) 

in very good agreement with 7 = 7 rcg = 43/32 = 1.34375. We have also tried to 
analyze our data assuming A = 1 . The quality of the fit is somewhat worse, although 
one could think that this is simply due to the additional neglected corrections to 
scaling^ that are still relevant at the values of iV we are working (the determination 
of A depends mainly on the data with N tot = 2000, i.e. on walks with N < 1000). 

We do not know how reliable the estimate ( |1.6|) is and a serious analysis of the 
systematic errors is practically impossible. In any case, without explicit assumptions 

1 An additional hint to A = 11/16 is the numerical observation of Barkema and Flesia j| that the 
average number of loops of length I forming a parallel contact scales as (m p )i ~ j-i.65io.Q5_ p rom 
this they estimate that the total number of parallel contacts behaves as (m p ) = a — b/N - 65 . On the 
other hand, (m p ) is proportional to cLzn(T) / dT if one includes an orientation dependent interaction, 
whence there should be generically a term ~ 7\r-o.65 w jy— n/ie m ZN (T). 

inclusion of a term with A = 11/16 worsens the quality of the extrapolation [ p2[ . 

3 A similar phenomenon occurs for SAWs on the square lattice jL8| . If one analyzes the end-to-end 
distance for short SAWs, with a single correction term, one finds A w 0.80. Only the inclusion of 
longer walks and more correction terms with Ai > 1 gives A sa 1.0. 
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on A, we can still obtain the lower bound 



7 > 1.3425 ± 0.0003. (1.7) 



Although this result is lower than the estimate ( |1.6| ), it clearly supports the fact that 
7 — 7re g - The prediction of Ref. |Tj|, 7 = 1.3385 ± 0.003 is instead clearly excluded. 



We should point out that if we analyze only data with small values of N, we would 
obtain a lower estimate for 7, in close agreement with the result (|1.4|) . This shows 
the crucial role played in this problem by corrections to scaling and the importance 
of performing simulations for very large values of N. 



2 The join-and-cut algorithm on the Manhattan 
lattice 

2.1 Description of the algorithm 

In this Section we will define the pivot and the join-and-cut algorithm on a Manhattan 
lattice. The Manhattan lattice is a two-dimensional square lattice on which bonds 
are directed in such a way that adjacent rows (columns) have antiparallel directions, 
corresponding to the traffic pattern in Manhattan. Although bonds are directed, there 
is no overall directional bias. Explicitly we will assume the following orientations: a 
vertical bond connecting the points of coordinates (x, y), (x, y+1) is directed upward if 
x is even, downward if x is odd; a horizontal bond connecting the points of coordinates 
(x, y) and (x + 1, y) is directed to the left if y is even, to the right if y is odd. 

Let us now define the pivot algorithm [23-^5]. It works in the ensemble of fixed- 



length walks with free endpoints: we will be interested in self-avoiding walks, but the 



algorithm can also be applied in a very efficient way to the Domb- Joyce model | 2q,p7 
to power- law walks [EEL and to interacting polymers far from the 0-transition [|2"E|,[ZU|. 

The algorithm works as follows [p5| . Given an iV-step SAW uj starting at the 
origin and ending anywhere, uj = {uj(0), . . . ,u(N)}, an iteration of the algorithm 
consists of the following steps: 

(i) choose randomly, with uniform probability, an integer k G {0, 1, . . . , AT — 1}; 

(ii) choose with probability P(g) an element g of the lattice point symmetry group. 
The probability must satisfy P(g) = P(g~ l ) to ensure detailed balance; 

(iii) propose a pivot move uj —> uj' defined by 

'(■\ _ / w 00 for 1 < z < fc, , 

U[i> \u(k)+g(u(i)-iu{k)) for k + l<i<N. 1 > 

The proposed move is accepted if uj' is self-avoiding; otherwise it is rejected and 
we stay with uj. 
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The probability P(g) must be such to ensure ergodicity: as discussed in Refs. [p5| , |30 
not all symmetry transformations are needed. In particular |30| the pivot algorithm 



is ergodic on a square lattice if P(g) is non-vanishing only for diagonal reflections, 
that is for reflections through lines of slope ±1. 

For the Manhattan lattice the point symmetry group is much smaller than the 
symmetry group of the square lattice. Indeed the lattice is symmetric only with 
respect to lines of slope —1 going through lattice points (x, y) with x + y even, and 
with respect to lines of slope +1 going through (x,y) with x + y odd. Therefore we 
will modify step (ii) in the following way: 

(ii) if ui x {k) +u) y {k) is even (resp. odd), let g be the reflection with respect to a line 
of slope —1 (resp. +1), going through uj{k). 

Since g 2 = 1, detailed balance is automatically satisfied. The tricky point is ergodicity. 
It can be proved by slightly modifying the proof of Section 5 of Ref . [[Kj . Using the 
same definitions, the basic observation is the following: if uu(k) = {xk,yk) is a walk 
site belonging to a diagonal support line of slope —1 (resp. +1), then Xk + yk is even 
(resp. odd). Using this fact the proof works without any change. 



We want now to define the join-and-cut algorithm [|TJ]. The tricky point here is 
that in general, given two walks on the Manhattan lattice, their concatenation is a 
walk that does not respect the bond orientation of the lattice. More precisely, given 
two walks iO\ , u>2 of lengths N\ and N2 , the concatenated walk u = uj\ o uj 2 respects 
the orientation of the lattice only ilfj 

mod^Aq) - u 2x (0), 2) = 0, mod(o; l2/ (Aq) - u 2y (0), 2) = 0. (2.2) 

This reflects the fact that only translations of two steps in each direction are symme- 
tries of the lattice. A second consequence of the bond orientation of the Manhattan 
lattice is the following: consider a walk u and cut it into two parts u>i and u 2 such that 
u — U\ ouj 2 - In general there is no lattice translation T such the translated walk Tuj 2 
respects the bond orientation of the lattice and satisfies (Tu 2 )(0) = c^i(0) = a>(0). 
Indeed this happens only if 

mod^O) - lu 2x (0), 2) = 0, mod(cu l2/ (0) - u 2y {0), 2) = 0. (2.3) 

With these two observations in mind we define our ensemble in the following way: 
T Ntot consists of all pairs (lji,u 2 ) of SAWs, each walk starting either in (0,0) or in 
(1, 1) and ending anywhere, such that the total number of steps in the two walks is 
some fixed even number N to t- Moreover we require the lengths of u)\ and uo 2 to be 
even. Explicitly 

iVtot/2-1 

T Ntot = (J (^(0, 0) U S 2k {l, 1)) x (S Ntot . 2k (0, 0) U S Ntot „ 2k {l, 1)) (2.4) 

k=l 

where Sk(x,y) is the set of k-step walks starting from (x, y) and ending anywhere. 
Each pair in the ensemble is given equal weight: therefore the two walks are not 

4 The only exception to this rule is when ui2 is a rod. 
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interacting except for the constraint on the sum of their lengths. Notice that there 
is a one-to-one correspondence between ^(0,0) and 5^(1,1) so that the ensemble 
defined in Eq. ( |2.4j ) is equivalent to 

iVtot/2-l 

T Ntot = (J S 2k (0,0) x S Ntot . 2k (0,0) (2.5) 

k=l 

One sweep of the algorithm consists of the following steps: 

(i) Starting from a pair of walks (001,002), we update each of them independently 
using some ergodic fixed-length algorithm. We use the pivot algorithm we 
described above. 

(ii) With probability 1/2 we interchange uo\ and oo 2 , i.e. (ooi,oo 2 ) — > (oo 2 ,ooi). 

(iii) We first check the parity of the endpoints: if mod^i^A^) — lj 2x (0), 2) = 1, we 
stay with {u)\, u 2 ) . Otherwise we attempt a join-and-cut move. We choose with 
uniform probability k G {1, . . . , N tot /2 — 1}. Then we concatenate the two walks 
001 and u 2 forming a new (not necessarily self- avoiding) walk c<j conc = uo\ o u 2 ; 
then we cut u CO nc creating two new walks u\ and co' 2 of lengths 2k and N tot — 2k. 
If oo\ and oo' 2 are self-avoiding we keep them; otherwise the move is rejected 
and we stay with u)\ and oo 2 . 

It is easy to see that the full algorithm is ergodic. 

The algorithm defined on the Manhattan lattice works essentially in the same 
way as the standard one. The only important difference is that one must check the 
parities of coi(Ni) and ^(O) before attempting the join-and-cut move. This check is 
successfull in 50% of the cases and thus the algorithm we have defined above should 
be essentially equivalent to the standard algorithm in which one performs two pivot 



updates of each walk for every join-and-cut move (in the notation of Ref. []14]] it 
corresponds to the algorithm with n piv = 2). In principle it is easy to avoid this 50% 
rejection by modifying the third step of the algorithm in the following way: 

(iii) [improved join-and-cut move]. If mod(a;x x (iVi) — u 2x (0),2) = 1, let u; conc = 
oj\ o Roo 2 , where Ruj 2 is the walk reflected with respect to the line of slope +1 
going through oo 2 (0); if mod(ooi x (Ni) — u 2x (0), 2) = 0, let uj couc =Wiow 2 . Then 
proceed as before. 

This modification attempts twice the number of join-and-cut moves with respect to 
the previous one, and thus it should be more efficient. However the difference in per- 
formance is not expected to be large. Indeed from the analysis of the autocorrelation 
times of Ref. |H , we obtain^ that the ratio r in t(npi v = 2)/ri nt (n P i v = 1) is equal to 



5 Notice that in Ref. Q Ti n t(« P iv) is reported in units of two iterations, each iteration consisting of 
^piv pivot attempts and one join-and-cut attempt. In order to make a correct comparison we should 
multiply T in t(n p i v ) by ?ipi v . Thus what wc report 8-s Tmt^piv — 2) corresponds to 2Ti n t(^piv — 2 

) of 

Ref. [pl. 
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1.6,1.5,1.2 for N tot = 500,2000,8000 respectively. For the algorithm on the Man- 
hattan lattice we expect r int (non— impr)/r int (imp) to be very similar to the previous 
ratio: therefore for N tot = 8000 we expect that the autocorrelation time for the origi- 
nal algorithm to be only 20-30% higher than the corresponding time for the improved 
one. One should also keep in mind that one iteration of the improved algorithm is 
slower than one iteration of the unimproved one, since one is doing twice the number 
of join-and-cut attempts. Therefore we expect that for N tot = 8000 the improved 
algorithm is only 10-20% better than the original one. Since the implementation of 
the improved algorithm required major changes of our codes, we decided to work with 
the unimproved version we described above. 



2.2 Determination of 7 from join-and-cut data 

Let us now discuss how the critical exponent 7 can be estimated from the Monte 
Carlo data produced by the join-and-cut algorithm. 

Let us start by noticing that the random variable N\, the length of the first walk, 
has the distribution 

7(N \ if #1 is even, 
,f(JVi) = { Z ^ tot) (2.6) 
if Ni is odd 

for 1 < N% < N tot — 1; here Z(N tot ) is the obvious normalization factor and is the 
number of iV-step walks going from the origin to any lattice point whose asymptotic 
behaviour for large N is given by (|1.2|) . The idea is then to make inferences of 7 from 
the observed statistics of N\. Of course the problem is that ( |1.2| ) is an asymptotic 
formula valid only in the large- N regime. We will thus proceed in the following way: 
we will suppose that ( |L2| ) is valid for all iV > N min for many increasing values of iV min 
and correspondingly we will get estimates 7(^V"tot, iV min ); these quantities are effective 
exponents that depend on N min and that give correct estimates of 7 as N min and N tot 
go to infinity. 

The determination of 7 from the data is obtained using the maximum-likelihood 
method (see e.g. [pi]] ). We will present here only the results: for a detailed discussion 



we refer the reader to [14 



Given N min} consider the function (from now on we suppress the dependence on 
AW 

e K .w=(l * (2. 7) 

Jv mm\ J- I n nflinrnnnn \ / 



otherwise 
and let X be the random variable 



X = log[N 1 (N tot -N 1 )] . (2.8) 

Then define 

X ccns (iV min ) = { f e 6Nn ^ (2.9) 
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where the average ( • ) is taken in the ensemble Tjv tot sampled by the join-and-cut 
algorithm. The quantity defined in Eq. ( [2.9|) is estimated in the usual way from the 
Monte Carlo data obtaining in this way X^f(N min ). Then 7(iV min ) is computed by 
solving the equation 

*MC (Wmin) = [X] thA (N min ), (2.10) 

where, for every function of N±, we define 

\f( N )] ( N ) - gISr fm N r\N tot - Nl p-\ 

here the sum is extended over even values of iV\. The variance of 7(iV min ) is then 
given by 

Var 7(iV min ) = , (2.12) 

U^;^JthM^mm)r 

where [X;X] = [X 2 ] - [X] 2 . We must finally compute Var (X^ s (N min )). As this 
quantity is defined as the ratio of two mean values (see Eq. Q2.9|)) one must take into 
account the correlation between denominator and numerator. Here we have used the 
standard formula for the variance of a ratio (valid in the large-sample limit) 

Ms) 

2.3 Dynamic behaviour 

Let us first discuss the simulations using the join-and-cut algorithm. We have per- 
formed high-statistics runs at N tot = 500, 2000, 8000 and 32000. In Table g we report 
the number of iterations and the CPU time per iteration on an Alpha-Station 600 
Mod 5/266. The total join-and-cut runs took 1 year of CPU on this machine. 

In Table ^| we report the acceptance fraction of the pivot and of the join-and- 
cut moves, and, for two different values of N m [ n , the autocorrelation times for the 
observable 

Y(N min ) = —- ^ - (2.14) 

that, according to Eq. fl2.13|) , controls the errors on 7. Notice that Y(2) = Xj (X) — 1, 
so that r intj y (2) = r int ,x- 

To compute the autocorrelation times we used the recipe of Ref. |32], Appendix C] . 
Indeed the autocorrelation function has a very long tail with a very small amplitude, 
due to the fact that r exp ^> r int y. For this reason the self-consistent windowing 
method of Ref. ^B, Appendix C] does not work correctly and underestimates the 



autocorrelation time. Following Ref. p2| , we have computed r int y as 



Tlnt.y = T int,Y + T ta ii,y. (2-15) 

Here Ti nt) y is the autocorrelation time computed using the self-consistent windowing 
method of Ref. [^5], Appendix C] with window W = 15r inti y and 

TtaU,Y = PY (W)W hg i^y^-J (2-16) 
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where py{t) is the normalized autocorrelation function and f p - lv is the acceptance 
fraction of the pivot algorithm. Of course Eq. (|2.16| ) is a very rough ad hoc estimate of 
the contribution of the tail. For X, it amounts to approximately 35% for N tot = 32000, 
24% for JVtot = 8000, 15% for N tot = 2000, and 5% for N tot = 500. Based on our 
experience with the pivot algorithm we have assigned to Tt a u,Y an error of 10%. The 
error on T- in t,Y was instead computed as in Ref. 



25 



We have first analyzed the acceptance fractions. The general analysis of Ref. 
predicts / piv ~ N~ p and f] c ~ N~ q , with p « 0.19 and q « 7 — 1 « 0.34. 
least-square fit to the data of Table || gives 



m 



A 



p = 0.19605 ±0.00001 
q = 0.32522 ±0.00004 



xl 
xl 



603, 
9197. 



(2.17) 
(2.18) 



The very large value of \ 2 indicates that the quoted errors on p and q, that are of 
purely statistical nature, are unreliable. Indeed a simple power law does not fit the 
data at this level of precision. A more realistic estimate of the errors would be 



p = 0.1960 ±0.0002, q = 0.325 ±0.004. 



(2.19) 



We have then performed an analogous analysis for the autocorrelation time T inti y(7v min ). 
A least-square fit of the form 

Tint,* ~ N z (2.20) 



gives 



0.847 ±0.006 
0.672 ±0.004 
0.747 ±0.003 



xl 
xl 
xl 



71 


for A min 


= 1, 


280 


for A min 


= 100 


50 


for A min 


= 200 



(2.21) 



Again, the very large value of \ 2 indicates that the error bars are underestimated by 
at least a factor of ten. The results are somewhat higher than the estimate reported in 
Ref. [14], z = 0.70±0.03. This is evident if one directly compares our results of Table 
|2] with the values of r int reported in Ref. [T4J (see footnote |): r int = 4.38(6), 9.47(21), 
20.0(6) for A" tot = 500, 2000, 8000 respectively. For A tot = 500 the autocorrelation 
times are similar, while for A tot = 8000 there is a factor-of-two difference. It is very 
difficult to believe that the join-and-cut algorithm has a different dynamic behaviour 
on the square lattice and on the Manhattan lattice. We have thus tried to understand 
if the authors of Ref. underestimated the autocorrelation times. Therefore we 
have determined T; nt from our data using the self-consistent windowing procedure 

4.166(2), 10.093(6), 



in 



with a window of 5ri nt> x as done in Ref 
29.13(2), 98.0(2) for N tot = 500, 2000, 8000 
gives z = 0.7271(2). This is now in good agreement with the results of Ref. 



We obtain r int 
32000 respectively. A fit of these results 



In conclusion we believe a fair estimate would be 



0.75 ±0.05 



(2.22) 



As noticed in Ref. |]I4" |, with a clever implementation of the algorithm, it is possible 
that the CPU-time per iteration increases as N a with o < 1. In particular it was 



10 



predicted that a = 1 — p where p is the exponent that controls the pivot acceptance 
fraction. In two dimensions we expect p = 0.1953(21) |25| and this is confirmed by 
our data on the acceptance fraction. From the data of Table |2] we obtain 

a = 0.865(1) xl = 55 (2.23) 

Keeping into account that the error bars are underestimated, we find reasonable 
agreement with the prediction a = 1 — p. Of course what one is really interested in 
is the autocorrelation times expressed in CPU units. 
We find 

Tint ~ N z+a ~ N 1 - 6 . (2.24) 

The algorithm is not optimal, but still represents an improvement with respect to 
other algorithms that behave as N~ 2 . 



3 PERM 

3.1 Description of the algorithm 

The pruned-enriched Rosenbluth method (PERM) is a chain-growth algorithm based 



on the well-known Rosenbluth (inverse restricted sampling) |33[ algorithm. In the 
latter, chains are grown step by step. In unbiased sampling, steps are randomly 
generated, and, if a new step violates the self- avoidance constraint, the configuration 
is discarded and one restarts from scratch. This leads to exponential 'attrition', i.e. 
the number of attempts needed to generate a chain of length n increases exponentially 
with n. 

In Rosenbluth sampling, such steps are replaced (if possible) by 'correct' steps. 
This diminishes the attrition problem, but it leads to a bias that has to be compen- 
sated by a weight factor. If m n is the number of free sites where to place the n-th 
monomer (i.e., the number of allowed steps at time n), then the weight of a chain of 
length iV is Wn oc II^li m n- The proportionality factor depends on the ensemble one 
wants to simulate, e.g. it is constant for the canonical ensemble with fixed fugacity. 
The main problem with the Rosenbluth method is the fact that these weights show 
large fluctuations for large N. Thus even large samples can be completely dominated 
by just a few events, which leads to large statistical errors. Even worse, for very large 
N, the part of the distribution of weights which should dominate might not have 
been sampled at all, and this can lead to systematic errors. More precisely, while 
the method is guaranteed to give an unbiased estimate of the partition sum, due to 
these fluctuations and to the convexity of the logarithm, estimates of free energies 
are systematically too small. 



The main idea in PERM |To| is to watch the weight Wn as N is growing, and to 
compare it with the estimated average weight for this value of N. If Wn is judged 
too big, a clone of the present configuration is made, both clones are attributed a 
weight Wn/2, and both are independently continued ('enrichment'). On the other 
hand, if Wn is too small, a (pseudo-) random number r is drawn uniformly from 
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[0, 1]. If r < 1/2, the configuration is abandoned ('pruning'), otherwise it is kept and 
its weight is doubled. This is most easily implemented by recursive function calls. A 
pseudocode is given in [|1^ . 

In this way the problem of large weight fluctuations is avoided, although not 
entirely. The set of generated configurations now is correlated due to cloning. Let 
us call a tour a set of configurations obtained by cloning from the same 'ancestor'. 
Although the weights of individual configurations vary between narrow limits, the 
weights of tours can fluctuate very much. If so, we are back at basically the same 
problem. This can be avoided by making additional biases. For instance, before each 
step one might look ahead k steps, and choose the direction accordingly ||34j| . This is 
rather time consuming for large k, nevertheless it can be efficient for dense (collapsed) 
systems. 

For diluted systems, as it is the case here, the efficiency can be improved using 
Markovian anticipation JT^,PH| . Here one keeps the last k steps in memory, and biases 



the next step on the basis of the statistics of (k + l)-step configurations accumulated 
during a previous run. 

On the Manhattan lattice, an iV-step walk with fixed starting point and fixed 
direction of the first step can be encoded as a binary sequence of length N — 1. A 
straight step is encoded as '0', while a bend by ±90° is encoded as '1'. Notice that one 
does not have to specify the direction of the bend, as only one bend is allowed at any 
time step. If the starting point has even x and y and the first step is upward, then the 
first bend is to the right (left) if it occurs at even (odd) times. After that, each bend 
is in the same (opposite) direction as the previous one if the number of in-between 
straight steps is even (odd). Let S = (s_fc, . . . , s ) = (s, s ) denote a binary string of 
length k + 1, and let Cjv,m(S) be the cylinder set of all N-step walks starting upward 
from (x, y) = (0, 0), for which steps m — k, . . . ,m are given by S. The total weight of 
this set in an unbiased ensemble is denoted by Wjv,m(S). Finally, we consider 

W N , m (S) 

PN,m(s \S) = = tTt 7~ 7T ■ I 3 - 1 ) 

1-4=0,1 W N,m{S, S'o) 

If SAW's were a k-th order Markov process, we would obtain perfect importance 
sampling for N-step walks if we would select the m-th step according to PN,m(so\ s )- 
Given the fact that SAW's are not Markovian, Eq. (|3.1|) comes closest to importance 
sampling given a finite (fc-step) memory. 

In practice, Eq. ( |3.1|) is not applicable as it stands, since it is practically impossible 
to acquire and store the needed statistics for all m and N. First of all we use the fact 
that PN,m(so\s) becomes independent of iV and m in the limit m, N — m ^> k. We 
therefore estimate it for large values of m and N by accumulating statistics for all 
N larger than some lower threshold (we used N > 400), and N — m fixed (we used 
N — m = 200). The fact that this p(sq|s) does not give correct importance sampling 
for m ~ N is not serious and does not reduce much the efficiency More serious is 
the fact that it would be inappropriate for small N and m. In particular, for m < k 
we have no string s to condition upon. We dealt with this problem by tapering the 
symbol sequence with a string of k zeroes, and by taking step so in the m-th step 
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with probability 

/ ,x const/m + lim mjV _ TO ^ 00 p 7V ' m (so|s) 
P m {so\s) = ^ f 77 -p 7 — rrr • (3-2) 

L S (,=0,1 [Const 771 + lim m ,JV-m— ooPJV,m(S0 sj 



with const « 20. 

Notice that all this fiddling is not crucial. The method per se is correct even if the 
probabilities for taking the m-th step were chosen badly, provided the probability is 
not set equal to zero for any allowed step. But efficiency can increase substantially 
with a good choice. Efficiency can be measured in several ways. The most straightfor- 
ward method consists in measuring tour-to-tour fluctuations. Since subsequent tours 
are statistically strictly uncorrelated, this gives immediately estimates for the error 
bars of any measured observable. In practice, not to waste CPU time, we bunched 
tours in groups of typically 10 4 tours, and measured average values and fluctuations 
over these bunches. 

A less direct but more instructive measure of efficiency is the number of tours in 
which a length N is reached in at least one configuration, normalized to the number of 
all tours. For other modified chain-growth algorithms such as incomplete enumeration 



or the Berretti-Sokal (BS) algorithm |37 l, it seems that this number decreases 



as const/ N, provided the cloning and pruning parameters are adjusted so that the 
total number of generated walks is independent of N. Improving the details of the 
algorithm cannot change this scaling law, but can change the constant considerably. 
For incomplete enumeration and the BS algorithm, the constant is 0(1). For PERM 
without Markovian anticipation const « 10. For Markovian anticipation with k = 19 
we found const ~ 130, see fig. [I]. 

A last check for efficiency is the following. As we said, the weights of tours can vary 
substantially, to the point that the distribution of these weights might be not properly 
sampled in the region which should dominate statistical averages. If this happens, 
the estimates of the free energies are too low, and all the errors are systematically 
underestimated. To show that this is not the case in the present simulation, we show 
in fig. ^| the measured distribution of tour weights for iV = 32000, together with the 
weighted distribution. We see that the maximum of the weighted distribution occurs 
at a weight which is still well within our sampling distribution. 

3.2 Determination of 7 using PERM data 

The raw data produced by PERM are estimates of the partition sum zn, for all 
N < N max . It is then straightforward to obtain estimates of 7. Assuming Eq. (|1.2| ). 
the exponent 7 can be estimated from ratios of z^. for three different values of iVj. 
Indeed, if one chooses two positive constants a < 1 < b, and defines 

X log Z a j\[ -\- y log Zfojy fog Z]\[ 

7eff N = ■ — - — , 3.3 

x log a + y log 



one checks easily that 7 e fr(iV) = 7, provided that 

x + y = 1 , ax + by = 1 . (3.4) 
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Figure 1: Probabilities to reach length N at least once in a tour, plotted against 
N. The curves are for PERM with /c-step Markovian anticipation, with k = 19, 10, 5, 
and 0. The straight line is 130/N. 




Figure 2: Full line: histogram P(\ogW) of tours with fixed weight W, on a logarith- 
mic horizontal scale. Normalization is arbitrary. Broken line: W x P(\ogW), again 
with arbitrary normalization. 
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If Eq. (|1.2|) is not fulfilled exactly, Eqs. ( |3.3j ) and ( j3.4|) define effective Af-dependent 
exponents. These exponents still depend on a and b. It is clear that statistical errors 
will be small if a <C 1 <C b, but then a wide range of chain lengths are needed. 
Since iV max is fixed and cannot be exceeded, this means that aN will be small, and 
systematic errors will become large if we want to reduce statistical errors. In addition 
to the ratio b/a, one has the product a x b at one's disposal. For fixed b/a, its optimal 
value depends on the way how the statistical error on increases with N. In the 
present case, they increase roughly as y/N. For practical reasons one also wants 
a -1 and b to be simple numbers. We found that good results were obtained with 
a = 1/2, b = 4 (and therefore x = 6/7, y = 1/7). 

The computation of the error bars is not trivial. Indeed one should take into 
account the fact that different zn are correlated and therefore compute the full co- 
variance matrix, which is practically impossible. Therefore we proceeded differently 
to obtain error estimates. For each bunch of tours we estimated [7 e ff(-^)]bunch accord- 
ing to Eq. ( fOl) , and from these values we obtained the variance of the latter. Notice 
that we did not use this procedure to obtain 7 e ff (iV) itself by taking the average over 
all bunches, as this would introduce a bias. 

In addition, in order to improve the estimates of X ccns (iV min ), cf. Eq. (f2.9p, ob- 
tained using the join-and-cut algorithm, we evaluated this quantity using the PERM 
estimates of z^. The estimate of the mean value is trivial. For the error bars, we 
evaluated [X cens (A^ min )]b U nch over bunches of tours and then estimated their variance. 
We used the same trick as for 7 e fj. Again, notice that we did not use this procedure to 
estimate the mean value itself, since this introduces a bias. We have checked however 
that this bias is extremely small. 



3.3 Dynamic results 

Let us now discuss the simulations using PERM. We have performed high-statistics 
runs using Markovian anticipation with k = 19. The dynamic data of the runs are 
reported in Table |3|. The total PERM simulation would have taken approximately 7 
months of CPU time on an Alpha-Station 600 Mod 5/266. 

We have first of all performed an analysis of the CPU time per independent 
configuration, which we expect to scale as aN z with z ~ 2. Fitting the data of Table 
0, and including in each fit only the data with A" > A" min in order to detect systematic 
effects, we have 

1.92 N min = 500, 

2.08 A" min = 2000, 1 J 

0.00003 A" min = 500, , , 

0.00001 A" min = 2000. [ ' ' 

The expected exponent z ~ 2 is clearly recovered. Let us now compare these results 
with the join-and-cut estimates. For the latter algorithm, the CPU time per inde- 
pendent configuration scales as (cf. tables [I] and ~ 0.0003A^ 16 and therefore the 
join-and-cut is faster in generating one independent configuration than the PERM as 
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long as iV > 1000. This is may not be however a fair comparison of the two algo- 
rithms. Indeed what we call an "independent configuration" for PERM is a bunch 
of walks which are correlated, but which still contain more information than a single 
walk. In order to make a fair comparison it is better to consider an observable and 
compare the CPU time needed to obtain the same error bars. We have therefore com- 
pared the results for our preferred observable X cens (iV m i n ), see table [|. We find that 
for N tot = 32000 both algorithms require essentially the same CPU time to produce 
data with the same error bars. Only for longer walks the join-and-cut algorithm is 
more efficient. Notice that this result is true if we use the Markovian anticipation. 
In the absence of this improvement, PERM would be about ten times slower and 
therefore less efficient than the join-and-cut algorithm already for iV pa 500. 



4 Estimate of 7 

Let us now discuss the results for the critical exponent 7. We will first analyze the 



data using the method presented in section |2.2| . This will allow us to use at the same 
time the data produced using the join-and-cut and the PERM. In table |] we report, 
for various iV m i n , the estimates of X^f(N min ) obtained using the two algorithms. 
We have checked that the raw data agree within error bars, therefore checking the 
correctness of our programs and error bars. 



In table |5] we report the estimates of the effective exponents 7 defined in Sec. 



A graph of these estimates together with the results of the fits described below is 
reported in fig. It is evident that the corrections to scaling are particularly strong 
and indeed 7 clearly increases with iV m i n and N tot . Under the only assumption that 
in the interval 1000 ^ N < 32000 the corrections to scaling have the asymptotic sign, 
i.e. if we exclude that 7 will decrease for larger values of N min and N tot , we obtain 
(using N tot = 32000, N mhl = 2000) 

7 > 1.3423 ±0.0018. (4.1) 

This clearly excludes the result of Ref. |]13|] . Under this very weak assumption, we can 
conclude that there is no evidence from our data of the non-universality predicted in 
Ref. 0. More conservatively, our data indicate that, if the effect is really there, it is 
extremely small: 

7re g -7 < 0.0014 ±0.0018. (4.2) 

We can improve this bound making additional assumptions on the corrections to 
scaling. For each value of iV m i n we assume that is well approximated by 

z N = Afi N N^ (l + aN- A ) (4.3) 

for N > N min . We are not able to determine A reliably and therefore we have 
performed fits for various fixed values of A. Our analysis works as follows. For each 
triple (iV m i n , iVtot, 0) and for each value of A, we can define an effective exponent 
7(-/Vmm, N to t, a) (from now on we suppress the dependence on A) using the natural 
generalization of Eqs. (|2.10 ) and ( [2.11 ) with a corresponding error A7(iV min , N tot , a). 
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Then for each iV m i n and A we determine the optimal value a op t of a and an estimate 
of the exponent 7(A^ min ), by minimizing with respect to a and 7(A^" min ) 

,1 2 



x 2 = E 

Aftot 



7(^V" m in, iVtot; 0) ~ l(N n 

A^(N min ,N tot ,a) 



(4.4) 



where the sum runs over N to t = 2000, 8000, 32000. The statistical errors are obtained 
in a standard fashion. 

In table ^| we report the results of the fits for various values of iV" min and for 
two values of A, A = 1 and A = 11/16. First let us observe that the estimate of 7 
decreases with increasing A. Therefore a lower bound independent of any assumption 
on the value of A is obtained using the results for A = 1. From table |6] we can estimate 

7 > 1.3425 ±0.0003. (4.5) 

We can try to do more and see if we can obtain from the data a rough indication 
of the value of A. In principle we can consider the results for various N m [ n and see 
for which value of A the estimates do not depend on A^ min . Notice however that the 
results with different iV min are strongly correlated — they are obtained essentially 
from the same data — and therefore an observed stability or upward trend should be 
interpreted with caution since it could simply be an effect of the correlations. If we 
consider the results of table ^|we see that the analysis with A = 11/16 is extremely 
stable: a opt and 7 are essentially constant and we can estimate 

7 = 1.3436 ±0.0003. (4.6) 

This result is in perfect agreement with the universal value 7 reg = 43/32. On the 
other hand the results for A = 1 show an upward systematic trend. However the 
effect is barely statistically significant (the results for iV min = 100 and N min = 500 
differ approximately by 1.5 combined error bars) and one could just think that the 
systematic increase is simply an effect due to the neglected corrections to scaling 
and/or a result of the correlations. Without further hypotheses, we cannot confidently 
go beyond the lower bound ( |4.5|) . However if we assume 7 = 7 reg , then A = 1 becomes 
less plausible and we can say that our data favour the presence of a non-analytic 
exponent. We cannot give a reliable estimate of A, but Saleur's value 11/16 fits our 
data very well. 

This conclusion is fully supported by the alternative analysis of the PERM data 
using Eq. (|3.3p. In fig. || we plot 7^, obtained with a = 1/2, b = 4, against jV~ n / 16 . 
If A has Saleur's value, and if there are no other corrections to scaling, we should 
expect a straight line intersecting the y-axis at 7. The most dramatic deviations from 
a straight line are strong period- four oscillations, observed also in |H| . Similar period- 



four oscillations are observed for SAWs on the square lattice |§ . They correspond to 
a singularity of the grand canonical partition sum at x — — . In fig. [| one may 

also observe a slight curvature which might suggest that A < 11/16. However, a more 
careful analysis, using also different pairs (a, b) and additional correction-to-scaling 
terms, suggests that this effect is not significant. In contrast, the fact that A < 1 
seems significant. Accepting for A a value between 1/2 and 0.7, we find again perfect 
agreement with 7 reg = 43/32, while the estimate of [O seems ruled out. 
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Figure 3: Estimates of the effective exponents 7 for N to t = 2000 (empty circle), 8000 
(empty square), 32000 (empty triangle) and estimates obtained with the optimal 
amplitude method for A = 11/16 (filled triangle) and A = 1 (filled circle). The 
dotted lines represent 7 reg = 43/32 and the estimate of Ref. |T3L 7 = 1.3385. 
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Figure 4: Effective exponents 7 e fj defined by Eq. ( |3.3|) , plotted against jV -11 / 16 . The 
extrapolation to the y-axis gives the estimate of 7. Only data with high statistics, 
i.e. iV < 1000 are included. 
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N tot 


iViter 


CPU-time 


500 


10 9 


0.248(1) 


2000 


12 • 10 8 


0.846(8) 


8000 


33 • 10 8 


2.41(3) 


32000 


23 • 10 8 


9.08(2) 



Table 1: Number of iterations and CPU times per iteration for various values of iV tot . 
CPU-times are expressed in ms and refer to an Alpha-Station 600 Mod 5/266. 



iVtot 


/piv /jc 


Tint,y(l) Ti n t,Y(100) 


500 


0.451943(11) 0.155468(11) 


4.912 ±0.031 8.171 ±0.023 


2000 


0.3446316(86) 0.0998629(77) 


14.30 ±0.23 19.38 ±0.12 


8000 


0.2624050(83) 0.0629925(65) 


49.6 ± 1.2 55.8 ±1.0 


32000 


0.199693(46) 0.039330(31) 


202.4 ± 7.4 206.7 ±7.0 



Table 2: Acceptance fraction for the pivot move (/ p i v ), for the join-and-cut move (/j C ) 
and autocorrelation times for the various values of N tot . Autocorrelation times are 
expressed in units of two iterations. 



N 

1 * max 


A'tour 


A^config 


A^indep 


CPU time 


500 


1.0 • 10 6 


1.140 • 10 6 


2.637- 10 5 


4.985 


2000 


8.876 • 10 8 


1.010 • 10 9 


6.533 • 10 7 


70.53 


4000 


1.756 • 10 8 


2.016- 10 8 


6.359 • 10 6 


298.6 


8000 


3.421 • 10 8 


4.035 • 10 8 


6.296 • 10 6 


1272.6(4) 


32000 


2.792 • 10 7 


3.703 • 10 7 


1.154 • 10 5 


244(7) • 10 2 



Table 3: Number of PERM tours (column 2) with N < A^ max , total number of 
configurations with A" = N max (column 3), number of independent configurations 
with A" = A^ max (i.e., of tours which reached A" = A^ max ; column 4), and CPU time 
per independent configuration (expressed in ms, on an Alpha-Station 600 Mod 5/266; 
column 5). 
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N tot = 


500 




N tot = 


2000 


■^v min 


yccns 
^PERM 


x^ccns 

A .ic 


■*V min 


yccns 
^PERM 


yccns 
A .ic 


2 


10.59979820(1214) 


10.5997734(851) 


2 


13.3738441(259) 


13.374033(132) 


20 


10.68281651(924) 


10.6828003(692) 


50 


13.4285517(215) 


13.428669(115) 


40 


10.75716591(695) 


10.7571316(555) 


100 


13.4796982(176) 


13.479777(100) 


60 


10.81801510(514) 


10.8179770(447) 


150 


13.5246592(144) 


13.524727(87) 


80 


10.86850452(365) 


10.8684807(356) 


200 


13.5644369(118) 


13.564528(77) 


100 


10.91057682(248) 


10.9105787(280) 


250 


13.5998096(96) 


13.599867(67) 


120 


10.94553280(170) 


10.9455605(216) 


300 


13.6313610(77) 


13.631408(58) 


140 


10.97427940(116) 


10.9742889(162) 


350 


13.6595411(61) 


13.659569(51) 


160 


10.99746609(70) 


10.9974737(115) 


400 


13.6847010(48) 


13.684742(43) 


180 


11.01557064(39) 


11.0155716(77) 


450 


13.7071218(38) 


13.707140(37) 


200 


11.02893162(21) 


11.0289259(46) 


500 


13.7270279(29) 


13.727030(31) 




JVtot = 


8000 




= 


32000 


7V min 


"V^ccns 

r H/lxlVl 


j^ccns 


N m in 


"V^ccns 


j^ccns 


2 


16.148415(85) 


16.148470(199) 


2 


18.92202(76) 


18.92211(36) 


100 


16.175459(77) 


16.175516(184) 


100 


18.92796(75) 


18.92804(35) 


200 


16.203131(70) 


16.203156(171) 


200 


18.93483(73) 


18.93486(34) 


300 


16.229231(64) 


16.229317(160) 


300 


18.94185(71) 


18.94184(33) 


400 


16.253729(58) 


16.253832(150) 


400 


18.94889(70) 


18.94884(33) 


500 


16.276733(52) 


16.276882(140) 


500 


18.95588(68) 


18.95582(32) 


600 


16.298361(48) 


16.298492(131) 


600 


18.96280(66) 


18.96275(31) 


700 


16.318720(43) 


16.318815(123) 


700 


18.96964(65) 


18.96960(31) 


800 


16.337910(39) 


16.337938(114) 


800 


18.97639(63) 


18.97635(30) 


900 


16.356014(36) 


16.356073(106) 


900 


18.98303(62) 


18.98297(30) 


1000 


16.373112(32) 


16.373133(100) 


1000 


18.98958(60) 


18.98951(29) 


1100 


16.389269(29) 


16.389252(93) 


1100 


18.99603(59) 


18.99599(29) 


1200 


16.404539(26) 


16.404505(87) 


1200 


19.00237(58) 


19.00235(28) 


1300 


16.418976(24) 


16.418930(81) 


1300 


19.00863(56) 


19.00853(27) 


1500 


16.445524(19) 


16.445475(71) 


1500 


19.02085(53) 


19.02071(26) 








2000 


19.04978(47) 


19.04973(24) 








2500 


19.07657(42) 


19.07649(22) 








3000 


19.10144(37) 


19.10130(21) 



Table 4: Raw data for Xy[c- We report separately the results obtained with the 
PERM and join-and-cut algorithms. 
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7(500) 


N min 


7(2000) 


■*V min 


7(8000) 


A^min 


7(32000) 


2 


1.32950 (3) 


2 


1.33263 (7) 


2 


1.33773 (21) 


2 


1.34059 (87) 


20 


1.32584 (5) 


50 


1.33423 (9) 


100 


1.33900 (25) 


100 


1.34097 (92) 


40 


1.32776 (7) 


100 


1.33516 (12) 


200 


1.33950 (29) 


200 


1.34105 (97) 


60 


1.32880 (9) 


150 


1.33569 (14) 


300 


1.33987 (33) 


300 


1.34105 (101) 


80 


1.32955 (12) 


200 


1.33609 (17) 


400 


1.34015 (36) 


400 


1.34106 (105) 


100 


1.33011 (15) 


250 


1.33638 (19) 


500 


1.34041 (40) 


500 


1.34112 (109) 


120 


1.33064 (20) 


300 


1.33661 (22) 


600 


1.34058 (44) 


600 


1.34122 (114) 


140 


1.33120 (28) 


350 


1.33681 (26) 


700 


1.34068 (48) 


700 


1.34136 (118) 


160 


1.33139 (40) 


400 


1.33692 (30) 


800 


1.34070 (53) 


800 


1.34143 (123) 


180 


1.33051 (63) 


450 


1.33705 (35) 


900 


1.34080 (58) 


900 


1.34139 (127) 


200 


1.33235 (132) 


500 


1.33711 (42) 


1000 


1.34077 (63) 


1000 


1.34143 (131) 










1100 


1.34075 (68) 


1100 


1.34162 (136) 










1200 


1.34073 (74) 


1200 


1.34172 (139) 










1300 


1.34072 (80) 


1300 


1.34147 (143) 










1500 


1.34077 (94) 


1500 


1.34140 (153) 














2000 


1.34232 (178) 














2500 


1.34256 (205) 














3000 


1.34250 (236) 



Table 5: Estimates of 7(iV tot ) for various values of N min and N tot . 



N min 


A = 11/16 


A = 1 




7 a opt 


7 aopt 


2 

100 
200 
300 
400 
500 


1.34306 (6) 0.43 (1) 
1.34300 (10) 0.51 (4) 
1.34318 (14) 0.54 (4) 
1.34344 (18) 0.57 (5) 
1.34357 (22) 0.59 (5) 
1.34356 (27) 0.59 (6) 


1.34159 (6) 0.85 (2) 
1.34200 (10) 1.75 (13) 
1.34198 (14) 1.93 (24) 
1.34192 (18) 1.98 (36) 
1.34230 (22) 2.23 (50) 
1.34251 (27) 2.40 (59) 



Table 6: Extrapolated estimates of 7 for A = 11/16 and A = 1. 
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